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Monte Carlo simulations have been performed for aqueous charged colloidal suspensions as a func- 
tion of charge density (a) on the particles and salt concentration C s . We vary the charge density in 
our simulations over a range where a reentrant solid-liquid transition in suspensions of silica and 
polymer latex particles has been reported by Yamanaka et al. [Phys.. Rev. Lett., 80 5806 (1998)]. 
We show that at low ionic strengths a homogeneous liquid-like ordered suspension undergoes crys- 
tallization upon increasing a. Further increase in o resulted once again a disordered state which is 
in agreement with experimental observations. In addition to this reentrant order- disorder transition, 
we observe an inhomogeneous to homogeneous transition in our simulations when salt is added to 
the disordered inhomogeneous state. This inhomogeneous to homogeneous disordered transition is 
analogous to the solid-gas transition of atomic systems and has not yet been observed in charged 
colloids. The reported experimental observations on charged colloidal suspensions are discussed in 
the light of present simulation results. 



PACS numbers: PACS nos.: 61.20.Ja, 82.70.Dd 



I. INTRODUCTION 



Mono disperse charge stabilized colloidal suspensions 
are studied with great interest because they exhibit a 
variety of structural orderings analogous to that occur in 
atomic systems 0, [2) . Unlike in atomic systems, the range 
and strength of the effective interparticle interaction 
among the constituent particles are tunable over a wide 
range in colloidal system. Hence, different structural or- 
derings can be realized with ease in these systems at am- 
bient conditions . Though the most dominant interpar- 
ticle interaction among like-charged particles in a charge 
stabilized colloidal suspension is screened Coulomb repul- 
sion, experimental observations in bulk suspensions such 
as reentrant transition J|, ||], vapor-liquid transition (J, |), 
stable voids coexisting with ordered or disordered regions 
H 0, ||, ^| have served as evidences for the existence of 
a long-range attractive component in the effective inter- 
particle interaction of like-charged colloids [|l0). Purely 
a repulsive pair-potential given by Derjaguin-Landau- 
Verwey-Overbeek(DLVO) theory can not explain the 
inhomogeneous nature or the phase separation observed 
in like-charged colloidal suspensions. The DLVO theory 
is advanced by Sogami-Ise by considering the counterion 
mediated attraction and regarding the macroionic part as 
a one-component system [p. 2|. The effective pair-potential 
U(r) p2| obtained from the Gibb's free energy is found to 
have a long-range attractive term in addition to the usual 
screened Coulomb repulsive term given by DLVO theory. 
Tata et al's computer simulations using this U(r) could 
explain the above mentioned experimental observations 
satisfactorily § | g ^ g. 

Recently there have been several alternative theories 
jl7], [i9) to understand the observed phase separa- 
tion in charged colloids. These theories, which are also 
known as "volume term theories" |[7| [ll| have strong 



similarity in approach but differ in technical details and 
are the extensions of Debye -Huckel theory for asym- 
metric electrolytes. These theories consider the total 
Helmholtz free energy of the system that arises from 
Coulomb interaction between all species in the system. 
The outcome of all these theories is that pair-potential 
between colloidal particles can remain repulsive as a func- 
tion of interparticle distance r, yet a gas-liquid or gas- 
solid phase separation can occur due to volume terms 
which are responsible for the spinodal instability. Re- 
cently Schmitzj2(| has brought out the connection be- 
tween volume term theories and the long-range attrac- 
tive term that is present in Sogami-Ise potential U (r) . It 
is shown that Sogami-Ise theory includes volume terms 
and the attractive term in U(r) arises from the electro- 
static interactions that are present in volume term the- 
ories. Further, the criticism on Sogami-Ise theory by 
Overbeek |2l|| and Woodward's]^] has been countered 
by SmalleyP^I and Schmitz[p^|. Hence, we use U(r) as 
the model potential in our MC simulations. 

The objective of this paper is two fold. (1) Carry 
out simulations as a function of a for volume fractions 
and salt concentrations which are the same as that used 
in experiments of Yamanaka et al |2l| and look for the 
reentrant solid-liquid phase behaviour observed by them 
. The reentered disordered state observed at high values 
of a is known to be inhomogeneous but the nature of the 
underlying phase transition has not been identified. (2) 
To identify conditions under which the reentered inho- 
mogeneous disordered state will transform to a homoge- 
neous state. Towards this end we perform a detailed MC 
study over a wide range of salt concentrations as C s is 
known to alter the position R rn of the minimum as well 
as the depth U m of U(r). Our simulations as a func- 
tion of a show the reentrant order-disorder transition in 
agreement with experimental observations and we go on 
to predict a inhomogeneous to homogeneous transition 
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as a function of C s , that is yet to be observed in charged 
colloidal suspensions. 

The details of MC simulations are given in the Sec. 
II. The effect of charge density on the phase behaviour 
of suspensions with different volume fractions at low salt 
concentrations are discussed in Sec. Ill, where the sus- 
pensions exhibit reentrant solid-liquid transition. Section 
IV includes the results on the influence of salt concentra- 
tion on the reentered disordered state. Comparison of 
our simulation results with the results of reported exper- 
iments is presented in Section V. A brief summary with 
conclusions are given in section VI. 

II. DETAILS OF SIMULATION 

Monte Carlo (MC) simulations are carried out using 
Metropolis algorithm with periodic boundary conditions 
for a canonical ensemble (constant TV, V, T where N, V 
and T are respectively the number of particles, volume 
and temperature). Particles of diameter <i(=I20nm) are 
assumed to interact via a pair potential U(r) having the 
functional form 

U[r) = 2 __ ^ Kd ) exp(- K r) (1) 

where A = 2 + nd coth(/td/2) and the inverse Debye 
screening length k is given as 

n 2 = ATte 2 {n p Z + C s )/{ek B T). (2) 

Ze is effective charge on the particle (related to the 
surface charge density by a — Ze/ird 2 ), C s the salt 
concentration, T the temperature(298K), e the dielec- 
tric constant of water and ks is the Boltzmann con- 
stant. The position of the potential minimum R m is 
given as Rm = {A + [A(A + 4)] 1 / 2 }/2k and its depth 
by U m = U(R m ). Both R m and U m depend on a and 
C s - For the required volume fraction </> (= n p ird 3 /6), 
the length I of the MC cell is fixed from the relation 
I = N/rip. By performing simulation with different 
number of particles, N, we found in our earlier simula- 
tions [l4|, [l6| that the results are the same within the 
statistical error for N > 432. Hence N is fixed at 432. 

The results presented in this paper correspond to 
0=0.005 and 0.03. At these volume fractions the 
deionised suspensions are known to crystallize into a 
body centered cubic (bec) ordering J26|, hence we have 
chosen particles placed on a bec lattice as the initial con- 
figuration in all our simulations. We monitor the to- 
tal interaction energy Ut and the first peak height of 
the structure factor, S max for identifying thermal equi- 
librium unambiguously |27|, ^8). Most of the simula- 
tions away from the transition (e.g., freezing) took ap- 
proximately 8 x 10 5 configurations to reach equilibrium, 
while those close to the transition needed nearly 6 x 



10 6 configurations. A MC Step (MCS) is defined as N 
attempted moves during which, on an average, each par- 
ticle gets a chance to move. The step size to move the 
particles during MC evolution process is chosen in such 
a way that the trial acceptance ratio is always around 
50%. After reaching equilibrium pair-correlation func- 
tion g(r) 7 coordinate averaged pair-correlation function 
g c (r) and S max are calculated using procedures reported 
earlier [^8[ ^9|. g c {r) is obtained by averaging the coor- 
dinates of the particles over a sufficiently large number 
of configurations [p9[ . Since g c {r) is free from thermal 
broadening and it helps in identifying the crystal struc- 
ture unambiguously. Further, for suspensions which ex- 
hibit solid- like behaviour, the g c (r) shows sharper peaks 
as compared to the corresponding g(r). The mean square 
displacement < r 2 (t) > for a chosen particle is defined as 
< r 2 (m) > = <| rj(m + n) — rj(n) | 2 >, where rj(m) is 
the position of i th particle after m MCS and < ... > de- 
notes the configurational average over the initial configu- 
rations n. Cluster size and its distribution are calculated 
as described in our earlier work|l4|. The total fraction 
of particles participating in the clustering F c , defined as 
the ratio of the total number of particles that partici- 
pate in the clustering to the total number of particles TV, 
has also been obtained for suspensions which exhibited 
inhomogeneous nature (phase separation). We identify 
the nature of suspension whether it is homogeneous or 
inhomogeneous by calculating the ratio d s /do, where do 
(= V3/2[l/n p } 1 ^) is the average interparticle separation 
calculated from the particle concentration n p and assum- 
ing a bcc-like coordination in the liquid-like order and d s 
is the average interparticle distance estimated from the 
first peak position in g(r). For a homogeneous suspension 
d s /do=l and for an inhomogeneous suspension d s /do < 
1. 



III. EFFECT OF CHARGE DENSITY 

We have carried out MC simulations for two volume 
fractions 4> = 0.005 and <f> = 0.03 as a function of charge 
density on the particles keeping the salt concentration 
at 2fiM. These parameters are same as that used by Ya- 
manka et a/.pq] in their experiments. Fig. 1 shows the 
pair-correlation functions and projections of the corre- 
sponding MC cell onto a zy-plane for suspensions with 
0=0.005 and cr=0.1, 0.23 and 0.5 /iC/cm 2 . For low values 
of a (=0.1 /iC/cm 2 ), g{r) and g c {r) showed a decay as 
a function of r, the corresponding particle positions(Fig. 
1(B)) showed disorder and the ratio d s /do is found to 
be one (see table I). These observations suggest the sus- 
pension at this value of a is homogeneous and liquid-like 
ordered. On the other hand for suspension with cr= 0.27 
/iC/cm 2 we find g(r) does not decay with r and the po- 
sitions of sharp peaks in g c (r) clearly indicate that the 
ordering is bec. The projection of particles in the MC 
cell shows perfect ordering and the calculated d s /do is 
found to be one confirming the bec ordered crystalline 
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structure is homogeneous. At high values of a ( = 0.5 
/iC/cm 2 ) g(r) and g c {r) are found to decay a function 
of r once again . Further, we observe the first peak in 
g(r) shifted to smaller r implying d s /do < 1. These ob- 
servations suggest that suspension is inhomogeneous and 
disordered. The corresponding projected particles in the 
MC cell also confirm this and showed particle free regions 
(voids) coexisting with disordered dense region. Since the 
g c {r) shows sharper features than the corresponding g{r), 
the dense disordered region is identified to be solid-like 
(amorphous). 

In order to identify the values of a at which the ho- 
mogeneous liquid (HL) freezes into a homogeneous crys- 
talline (HC) order and then to an inhomogeneous dis- 
ordered (also referred as a phase separated (PS)) state, 
the structural parameter S ma x is calculated as a func- 
tion of a and is shown in Fig. 2. The sudden increase 
in S max at er=0.23 /iC/cm 2 (Fig. 2(a)) corresponds to 
the freezing transition and the sudden drop in S max at 
tr=0.33 /iC/cm 2 corresponds to the transition from ho- 
mogeneous crystalline state to a phase separated (gas- 
solid coexistence) state. Simulations for suspension with 
0=0.03 showed similar phase behaviour except that the 
corresponding transitions have been found to occur at 
lower values of a (see Table I). 

These observations are understood from the depen- 
dence of U (r) on a. At low charge density the suspensions 
remain homogeneous because the particles which are at 
a distance do will always experience screened Coulomb 
repulsive interaction as do < Rm (see curve a of Fig. 
3). The structural ordering in the homogeneous state 
depends upon the strength of U(r) at r = do- As the 
charge density on the particles is increased the strength 
of repulsive interaction increases. This results in freez- 
ing into a homogeneous crystalline state. Since R m de- 
creases monotonically with increase in a, particles which 
are at separation of do will experience attraction when 
R m < d (see curve 6). Due to strong attraction the 
particles condense into a solid-like dense phase leaving 
some particles in the rare phase (gas-like) which appear 
as voids when the fraction of the volume occupied by the 
rare phase is smaller than that of the dense phase. Thus 
the occurrence of freezing of a homogeneous liquid into 
a homogeneous crystalline state at lower values of a and 
the system reentering into a disordered inhomogeneous 
phase (gas-solid coexistence) at higher values of a are 
understandable. 



IV. EFFECT OF SALT ON INHOMOGENEOUS 
DISORDERED PHASE 

There have been reports of theoretical |29|, |3^] 
and experimental investigations^, |2(| [32| on the effect of 
C s on melting/freezing of a homogeneous colloidal sus- 
pension. However, the conditions under which an in- 
homogeneous disordered state, which occurs in highly 
charged colloids, becomes homogeneous are not known. 



When the well depth U m < ksT, the attraction between 
particles becomes weak, hence the suspension can remain 
homogeneous. It can be seen from Fig. 3 (curve d) that 
the well depth U m decreases to values smaller than ksT 
at high salt concentrations. So, we performed MC simu- 
lations over a wide range of salt concentrations keeping 
the other suspension parameters fixed and results are 
summarized in Table II. Fig. 4 shows the effect of in- 
creasing salt concentration on the disordered dense phase 
which coexisted with voids at C s — 2/iM (Fig. 1(c)). Fol- 
low values of C s , g{r) shows a decay, but the structural 
correlations persist to the full length of the MC cell (Fig. 
1(c)) suggesting that the disordered dense phase is con- 
nected. The fraction of volume occupied by the voids 
is found to be smaller than that occupied by the dense 
phase. Hence the voids constitute the minority phase 
at low values of C s . The particle positions in the MC 
cell(see Fig. 1(c)) confirm this. When the salt concen- 
tration is increased we observe structural correlations in 
g(r) (Fig. 4(a)) to extend up to a few times the diam- 
eter of the particle. The first peak height is also found 
to be very high and the first peak position matches ex- 
actly with R m . These observations imply the formation 
of dense phase clusters with R m as the average distance 
between particles within the cluster. The other sharp 
peaks appearing at distances beyond the first peak cor- 
responds to short-range ordering of particles within the 
clusters (intracluster ordering). 

Since the clustering occurs at large values of C s , the or- 
dering of particles within a cluster could be liquid-like or 
solid- like (amorphous). In order to distinguish between 
the two types of ordering we have calculated g c (r), g(r) 
and the mean square displacement < r 2 (t) > of particle 
within the cluster. We compare the behaviour of these 
quantities with that corresponding to a vapor-liquid co- 
existence state. This vapor-liquid coexistence state ob- 
tained by performing simulations for suspension param- 
eters at which a vapor-liquid coexistence has been re- 
ported |IL |J H . It can be seen from Fig. 5 that the peaks 
in g c (r) (curve a) corresponding to C s =550/iM are much 
sharper with increased peak height as compared to the 
corresponding g(r) (Fig. 4(a)). Whereas the g c {r) cor- 
responding to the vapor-liquid coexistence (dotted curve 
in Fig. 5) is much lower in peak height and also poorer 
in features as compared to the corresponding g(r) (inset 
of Fig. 5). The mean square displacement of particles 
inside the cluster corresponding to C s =550//M shows a 
saturation behaviour as a function of MC time and un- 
dergo much smaller displacement as compared to that of 
particle within a liquid-like ordered cluster. These obser- 
vations unambiguously suggest that the clusters formed 
at higher values of C s for highly charged colloids are 
solid-like. These clusters can be observed either in the 
scattering or microscopy experiments provided the den- 
sity of solvent is matched with the density of colloidal 
particles. Tata et al have observed the coexistence of 
voids with dense phase disordered regions ^ |j| and ex- 
istence of clusters under density matched conditions M 
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in highly charged polycholorostyrene sulfonate particles 
(PCSS) dispersed in an aqueous medium. Density match- 
ing for PCSS particles has been achieved by redispersing 
them in 60% glycerol aqueous solution. 

Simulations at a salt concentration of 1000 fxM showed 
only a single peak in g(r) (see Fig. 4(b)) at r = R m . The 
corresponding projection of particles in the MC cell (see 
inset of Fig. 4 (b)) do not show clustering and are dis- 
tributed uniformly. These observations suggest that the 
suspension with C s = 1000 /iM is homogeneous and non- 
interacting (gas- like). Thus we observe for C s < 550 /iM 
an inhomogeneous state in the form of dense phase clus- 
ters having solid-like ordering and a homogeneous gas- 
like disordered state at C s =1000/iM. Hence the inhomo- 
geneous (PS) to homogenous disordered transition is ex- 
pected to occur in between these salt concentrations. We 
identify this transition by monitoring the total fraction 
of particles participating in the clustering F c (Fig. 6(a)) 
and the first peak height g ma x of g(r) (Fig. 6 (b)). Note 
the sudden change in both of these parameters around 
C s = 570/iM as C s is increased. This change is associ- 
ated with the inhomogeneous (PS) to homogenous (HG) 
transition or in other words transition from a gas-solid co- 
existence to a gas phase. The large values of g m ax in the 
inhomogeneous state arises due to scaling of g(r) with re- 
spect to particle concentration n p . The dense phase con- 
centration rid is estimated from the first peak position ( = 
R m ) of g(r) and is found to be several times higher than 
n p (see Table II). In the case of suspension with 0=0.03 
the transition is found to occur at C s =270/iM which is 
lower than C s = 570/iM for tj> = 0.005. Suspensions of 
higher volume fraction undergo PS to HG transition at 
lower values of C s . This is due to increased screening 
of particles in suspensions of higher volume fraction as 
compared to that in suspensions with lower cj). 



V. COMPARISON WITH EXPERIMENTS 

It is important to compare our simulations results with 
that observed experimentally. Our simulations as a func- 
tion of a revealed the existence of the three phases viz., 
a homogenous liquid (HL) at low a, a homogenous crys- 
talline (HC) state at intermediate values of a and an 
inhomogeneous disordered state (PS) at high values of a. 
Further, we showed that the homogenous liquid freezes 
into a homogeneous crystalline state at about a > 0.22 
/iC/cm 2 and disorders once again at about a ge 0.33 
/iC/cm 2 . Yamanaka et aL[^H| have reported freezing of 
a homogenous liquid-like ordered aqueous suspension of 
charged particles with 0=0.005 into a homogenous crys- 
talline state for charge densities beyond 0.24 /iC/cm 2 . 
The suspension is found to disorder once again when a 
is increased beyond 0.40 /xC/cm 2 . This reentered disor- 
dered state is identified to be inhomogeneous from the 
scattering measurements. The measured average inter- 
particle separation d exp is found to be less than c?q |25j . [j3"| , 
which implies the nonspacefilling nature of the suspen- 



sion at high charge densities. Thus present simulations 
are in good agreement with the experimental observa- 
tions. Further, our simulations at higher volume fraction 
have revealed that the values of a at which HL to HC 
transition and HC to PS transition occur are lower as 
compared to that in suspensions of lower <j). This result 
is also in agreement with experimental observations ]25| . 

Present simulations on the reentrant inhomogeneous 
disordered state as a function of salt concentration re- 
vealed that the inhomogeneous state undergoes a tran- 
sition to homogeneous gas phase beyond a critical salt 
concentration. This gas-solid transition which can be 
observed in highly charged colloids by varying salt con- 
centration, has not yet investigated experimentally. Fur- 
ther, we notice from these simulations the absence of 
vapor-liquid coexistence for highly charged colloids. In 
this context, it is worth mentioning the result from MC 
simulations on a system of hard particles with an at- 
tractive Yukawa interaction by Hagen and Frankel 
They showed the disappearance of vapor-liquid coexis- 
tence when the range of attractive part of the Yukawa 
potential is less than approximately one sixth of hard- 
core diameter. In present simulations the range of at- 
tractive and repulsive parts of U (r) have been varied by 
changing a and C s . The range of attraction is found to 
be smaller for suspensions of high charge density parti- 
cles as compared to the suspensions of low charge density 
particles. Hence, one can observe a vapor-liquid coexis- 
tence in suspension of low charge density particles ||, [j| 
but not in suspensions of highly charged particles. 



VI. CONCLUSIONS 

MC simulations using a pair-potential with long-range 
attractive term show a homogenous liquid, homogenous 
crystal and a disordered inhomogeneous phase as the 
charge density is varied, suggesting a reentrant order- 
disorder transition, which is in agreement with the ob- 
served experimental results. The reentered inhomoge- 
neous phase is identified as the gas-solid coexistence state 
and occurs due to strong attraction experienced by the 
highly charged particles. Suspensions exhibit homoge- 
neous liquid-like and homogenous crystalline orders at 
lower values of a due to the dominant screened Coulomb 
repulsion at average interparticle separation. The inho- 
mogeneous disordered state undergoes a transition to a 
homogenous gas state upon increasing the salt concen- 
tration. This transition is analogous to the gas-solid 
transition observed in atomic and colloid-polymer mix- 
tures as a function of temperature and free polymer, 
respectively p5|. We believe that our prediction of gas- 
solid transition driven by the salt concentration would 
motivate experimentalists to investigate highly charged 
colloids in search of this transition as well as for detailed 
phase diagrams with a gas-liquid or gas-solid coexisting 
phases. 
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FIG. 1: (A) g(r) vs r for different charge densities at a salt 
concentration of 2/iM and <f> = 0.005. Curves a, b and c corre- 
spond to a = 0.1, 0.23, 0.50 /iC/cm 2 respectively. The inset 
shows (/ c (r) vs r. The curves a', b' , c' correspond to the same 
parameters as those of curves a, b and C. Curves 6, c, b' and 
c' are shifted vertically for the sake of clarity. (B) Projec- 
tion of time-averaged particle coordinates in the MC cell for 
parameters same as a, b and c. 

FIG. 2: Structural parameter S ma x as a function of a for a 
suspension with C s =2.0 /uM and <j> =0.005. The line drawn 
through the points is a guide to the eye. 

FIG. 3: Pair potential U(r)/kBT for different suspension pa- 
rameters with <^)=0.005. Curve a and b correspond to a= 0.15 
and 0.5 /u,C/cm 2 , respectively and C s =2.0 fjM. Curve C and 
d correspond to C s =550.0 and 1000 fjM, respectively and a 
= 0.15 /xC/cm 2 . The vertical line corresponds to the average 
interparticle separation do. 

FIG. 4: g(r) vs r at two different salt concentrations (a) 
C s = 550.0 ,uM and (b) at C s = 1000.0 fxM. Insets show the 
projection of particles in the MC cells corresponding to the 
two values of C 3 . Other parameters of the suspension are <j> 
=0.005, a =0.5 fiC/cm 2 . 

FIG. 5: (a) g c (r) vs r for suspensions exhibiting gas-solid 
(full curve: <f>— 0.005, a= 0.5 /iC/cm 2 ) and vapor- liquid 
(dotted line: <j>= 9.02 x 10~ 3 , a= 0.21/iC/cm 2 , C s =40.8mM, 
d=0.109/im) coexistence. g(r) vs r for the vapor-liquid co- 
existence is shown as inset, (b) Mean square displacement 
< r 2 (t) > vs t (in units of 1000 MCS) for particles within a 
solid-like (full curve) and liquid-like (dashed curve) clusters. 

FIG. 6: Quantities (a) F c and (b) g m ax as a function of C s 
used for the identification of the transition from a inhomo- 
geneous disordered state to homogenous gas phase. Other 
parameters of the suspension are same as that of Fig. 4. 
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TABLE I: Inverse Debye screening length (k), position of the 
potential minimum (R m ), depth of the potential well (U m ), 
average interparticle separation (do) and average interparticle 
separation d s calculated from first peak position of g(r) for 
different values of a corresponding to two suspensions with 
widely differing volume fractions. The abbreviations HL, HC 
and PS represent the homogeneous liquid, homogeneous crys- 
talline and phase separated states, respectively. 

4> a (/i C/cm 2 ) Kd Rm/d U m /kBT d /d d s /d State 



0.005 





.065 





.667 


7 


,350 


-0.088 


5. 


143 


1 


.0 


HL 


0.005 





.100 





.718 


6 


.841 


-0.221 


5. 


.143 


1 


.0 


HL 


0.005 





.150 





.786 


6 


,270 


-0.539 


5. 


143 


1 


.0 


HL 


0.005 





.200 





.848 


5, 


,837 


-1.023 


5 


143 


1. 


.0 


HL 


0.005 





.230 





.883 


5. 


,600 


-1.399 


5. 


.143 


1 


.0 


HC 


0.005 





.273 





.930 


5. 


,340 


-2.050 


5. 


143 


1 


.0 


HC 


0.005 





.300 





.960 


5, 


,190 


-2.547 


5 


143 


1. 


.0 


HC 


0.005 





.410 


1 


.070 


4, 


,690 


-5.173 


5. 


.143 





.9 


PS 


0.005 


0. 


.500 


1 


.152 


4. 


,383 


-8.119 


5. 


143 





.8 


PS 


0.030 





.060 


1 


.021 


4. 


,890 


-0.107 


2. 


.830 


1 


.0 


HL 


0.030 





.150 


1 


.462 


3. 


,542 


-0.850 


2, 


830 


1 


.0 


HC 


0.030 





.200 


1 


.660 


3 


,100 


-1.608 


2, 


830 


1 


.0 


HC 


0.030 





.273 


1 


.900 


2. 


,841 


-3.142 


2, 


830 


1. 


.0 


HC 


0.030 





.300 


1 


.990 


2. 


,744 


-3.853 


2, 


830 


1 


.0 


PS 


0.030 





.410 


2 


.300 


2 


,460 


-7.348 


2. 


830 





.9 


PS 


0.030 





.500 


2 


.530 


2 


,300 


-10.90 


2. 


830 





.8 


PS 



TABLE II: Inverse Debye screening length (k), position of the 
potential minimum (R m ), depth of the potential well (U m ), 
and the ratio of dense phase concentration rid to homoge- 
neous particle concentration n p for different values of C a for 
suspensions with volume fractions 0=0.005 and 0.03 exhibit- 
ing inhomogeneous and homogeneous states. The abbrevia- 
tions PS and HG represent the phase separated state and a 
homogeneous gas state, respectively. 

4> cr (/J, C/cm 2 ) C s {pM) nd R m /d Um/ksT nd/n p State 



u.uuo 


u.o 


i n 

1U 


1 RP17 

1 .UU 1 


9 97n 

O.Z 1 u 


-y .yuy 


3 8QD 
o.oyu 


PS 


n nn^ 


u.o 




9 0^7 


9 7nn 

z. / uu 


1 7fi 


fi Q8D 
u.you 


PS 


u.uuo 


u.o 


ou 


9 07^ 


9 D7D 
z.u / u 


i n c.n 

- 1U.OU 


10.00 


PS 


u.uuo 


u.o 


1UU 


4 n8fi 


1 74n 


8 ^41 

-0.04:1 


97 1 1 


PS 


n nn^ 


u.o 




^ fi88 
o.uoo 




^ 7^7 
-o. t o t 


D4 


PS 

1 o 


u.uuo 


u.o 


ouu 


fi Q^fl 
u.you 


1 499 


4 9R7 


4Q 71 
4y. / 1 


PS 


u.uuo 


n ^ 
u.o 


^n 

OOU 


7 474 


i.oyu 


^ 774 
-O. 1 1 4 


oo. iy 


PS 

-TO 


n nn^ 

u.uuo 


u.o 


4nn 

4UU 


7 Qsn 

J .you 


i ^fin 

1 .ouu 


3 ^89 
-o.ooz 


oo. yu 


pc 
1 o 


n nn^ 

u.uuo 


u.o 


ouu 


o.yuo 




9 SfT 1 ! 


fif) 79 
uu. t z 


PS 


n nn^ 

u.uuo 


u.o 


^n 

OOU 


Q ^4n 
y .oiu 


1 .OlO 


9 ^8^ 
-^.ooo 


fi^ ^4 

UO.04 




n nn^ 

u.uuo 


u.o 


Ron 

uuu 


Q 74"^ 

y. 1 4:0 


i ^nn 

1 .ouu 


9 


i nnn 

1 .uuu 


HP 

1 1 V 1 


u.uuo 


u.o 


1 uu 


1 fl ^91 


1 978 
i .z / o 


9 no^ 


i nnn 

1 .uuu 


HP 
1 1 V 1 


0.005 


0.5 


800 


11.241 


1.261 


-1.860 


1.000 


HO 

1 1 V.I 


0.005 


0.5 


900 


11.918 


1.246 


-1.670 


1.000 


HG 


0.005 


0.5 


1000 


12.558 


1.234 


-1.524 


1.000 


HG 


0.030 


0.3 


10 


2.285 


2.474 


-3.933 


1.500 


PS 


0.030 


0.3 


50 


3.389 


1.915 


-3.548 


3.230 


PS 


0.030 


0.3 


100 


4.395 


1.680 


-2.858 


4.782 


PS 


0.030 


0.3 


150 


5.211 


1.565 


-2.335 


5.916 


PS 


0.030 


0.3 


200 


5.915 


1.495 


-1.958 


6.786 


PS 


0.030 


0.3 


250 


6.544 


1.446 


-1.680 


7.500 


PS 


0.030 


0.3 


300 


7.117 


1.410 


-1.471 


1.000 


HG 


0.030 


0.3 


400 


8.145 


1.358 


-1.176 


1.000 


HG 


0.030 


0.3 


500 


9.055 


1.323 


-0.980 


1.000 


HG 


0.030 


0.3 


600 


9.883 


1.296 


-0.841 


1.000 


HG 



250 
200 
3 150 
off 100 
50 




0.08 








5 
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o 0.06 



V 



0.04 



0.02 



0.00 
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